Dynamics of an inchworm nano— walker 
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An inchworm processive mechanism is proposed to explain the motion of dimeric molecular 
motors such as kinesin.We present here preliminary results for this mechanism focusing on 
observables like mean velocity, coupling ratio and efficiency versus ATP concentration and 
the external load F. 

Molecular motor proteins transform the energy of ATP hydrolysis into mechanical work per- 
forming discrete steps along a periodic track. The experimental work on protein motors has 
stimulated a wide variety of modelizations, most of them based in ratchet-like potentials. The two 
main candidates for the walking mechanisms of dimeric motors were inchworm or hand-over-hand. 
In the first case it is assumed that the first leading head advances one step which is followed instan- 
taneously by the trailing head. In the hand-over-hand mechanism the second head advances two 



steps overpassing the first head. In Ref. |2|] some experimental evidence was presented which seem 
to support the inchworm mechanism. Nevertheless, more precise experiments show that myosin-V 
walk in a hand-over-hand way These two different mechanisms imply different conformational 
changes in the protein structure during ATP hydrolysis. Moreover they imply a different response 
with respect the experimental control parameters [ATP] and F. 

Although it is commonly accepted now that some processive members of kinesin, myosin and 
dynein families seem to walk in a hand-over-hand fashion, it is still worth analyzing the inchworm 
mechanism, which could hold for other type of motors. For these reasons, we will present here a 
very simple model walking in a inchworm fashion with parameter values in the biological scale. We 
will also focus on some implications with experimental relevance. 
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It was showed in Ref. 1] that kinesin uses a single ATP molecule to perform each step. Such relation 
is called the coupling ratio, which for low external loads is about 1. The temporal distribution of 
these steps is random due to the ATP diffusion until it reaches the motor. After binding the 
nucleotide, hydrolysis and the consequent conformational change take place displacing the whole 
motor a certain distance which is usually equal to the periodicity of the track. All this process can 
occur even in the presence of an opposing external force F and at low ATP concentration, although 
both regimes decrease the mean velocity. 

The inchworm walking mechanism can be modeled as two linearly-coupled particles interacting with 
a ratchet potential. The conformational cycle is introduced as a stretching and posterior relaxing 
of the coupling spring. This way of modeling was introduced in Ref. [J] showing that thermal 
fluctuations are not strictly necessary in order to achieve the motion. Other works P, |(], Q, 0, |l| 
are also based on this approach. The main difference between them is the way they model the 
mechanical changes under the input of chemical energy. While 4] considers the conformational 
change as an increase of the equilibrium length of the spring, other literature introduce asymmetric 
frictions or switches on the ratchet potential. Here we will explore a different and simple way which 
allows to control with precision the amount of the input energy. Furthermore, we apply a kinetic 
methodology 10] based on enzymatic inhibition to get an analytical expression for the velocity as a 
function of the ATP concentration and the external force F. Finally we analyze the coupling ratio 
and the efficiency at different values of F. 

We will consider the motor as two particles coupled by a spring. The set of equations in the 
overdamped limit are, 

Aii = -V'( Xl ) - k(x x -x 2 -L)- f s (t) - | + £i(t) 

Xx 2 = -V'{x 2 ) + k(x 1 - X2 -L) + f s (t) - | + &(t), (1) 

where x%,X2 are the position of the trailing and the leading head, respectively, k is the stiffness of 
the harmonic spring with equilibrium length L, which is also the periodicity of the ratchet potential 
V(x). Such potential has an asymmetric factor a and a barrier height Vq (See Fig^). f s (t) is the 
random chemical force and F is the external load. The thermal force is emulated through a zero 
mean Gaussian-white noise with a correlation, = 2^kBT5(t — t'). A is the friction and 

ksT is the thermal energy. Then, at thermal equilibrium the two particles will lay, most of the time, 
on two consecutive potential minima. We assume that at random intervals of time, an energetic 
nucleotide like ATP will bind the motor and a stretching force f s will act on the system until the 
total length of the motor will be doubled, i.e. x 2 — x\ = 2L. If E is the hydrolysis energy of the 
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FIG. 1: a) Scheme of the ratchet potential the positions of the leading head (black particle) and the trailing 
head (gray particle) at the three stages of the motor, i) Rest configuration, ii) The motor at the end of the 
stretching, iii) Final stage, b) Simulated mean velocities (v) of the center of mass versus p for F — 0, 1, 3 pN 
drawn with circles, squares and triangles, respectively. Dashed lines are Michaelis-Menten fits. The resulting 
kinetic parameters are shown in the table. 

molecule and L is the displacement that it is performed, then we take the chemical force f s = E/L. 
The values of the parameters have been chosen in a nano scale to mimic some molecular motors 
such the kinesin. The periodicity of the potential L is taken to be the periodicity of microtubules, 
8nm. The asymmetric factor a = 0.8 and Vq = bOpNnm optimizes the efficiency of our model. 
E = WOpNnm corresponds to an accepted value for the energy of hydrolysis of an ATP and 
thermal energy is ksT = A.lpNnm. The stiffness of the motor is chosen k = lpN/nm and the drag 
force A = 2 • 10 ~ pNs/nrn. When the motor is free from ATP, p € (0, 1) is the uniform probability 
per time step At of binding one molecule. When it occurs, more ATP binding is forbidden and 
stretching takes place until the elongation is 2L. Then, the stretching force disappears and the 
spring relaxes. When %2 — x% is again L, one cycle is completed and ATP binding is allowed. In 
the absence of external load, this mechano-chemical cycle induces a L displacement of the motor 
towards one end of the potential. Fig. shows the scheme of the process. On the other hand, 
our approach controls how much energy Ej- is applied to the system by simply multiplying E 
by the number n of ATP consumed: Ej- = nE. The mean velocity of the motor, when the ATP 
concentration is saturant and F = 0, is maximum and dependent only on the intrinsic properties 
of the motor and by E. Let t on be the time spent to perform a single step, i.e. the stretching plus 
the relaxing time. Thus, V ma x = L/t on . Using the given values of the parameters, simulations show 



that t on ~ 0.012s, which gives V max ~ 667nm/s. However, the global speed (v) will be slowed 
down when the ATP concentration decreases. Typically, the [ATP]-dependence on (v) is given 
by the Michaelis-Menten relation In our model, we have previously defined p as the uniform 
probability to get an ATP per time step At and with p = while the motor stretches and relaxes. 
It can be accepted that, as the reaction frequency is proportional to [ATP], and then [ATP] is 
proportional to p. Then, we have 

<"» = v "-"Kihr P ' (2) 

where Km is the Michaelis constant for the probability. From now on, we will deal with p and not 
with [ATP]. Fig^> shows how the michaelian behavior fits well the simulated values of the mean 
velocity. However, for finite values of P, both kinetic parameters V max and Km change. In Ref . [lol| 
it is shown that the effect of the external load in kinesin can be interpreted as an inhibition process. 
This introduces a P-dependence on the two kinetic parameters 

W^) = 7^ (F r Q) Km(F)=K m (F = 0) ] + gggg . ( 3 ) 

1 + K tu {F s /F-l) 1 + K iu (Fs/F-l) 

and allows to express the velocity of the motor as a function of the two control variables p and P. 
Ps is the stall force, i.e. the maximum load that the motor is able to carry. In our simulations, 
Fs ~ 5.25pN. Ki u and Ki c are the uncompetitive and competitive inhibition constants, respectively, 
and are a quantitative measure of how P affects the motor when it is free from nucleotide (Ki c ) or 
when it has an ATP (Ki U ). From table values of Fig^) we fit the values of the inhibition constants 
obtaining K iu ~ 0.338 • 10~ 6 and K ic ~ 2.131 • 10" 6 . As they are dissociation constants, the effect 
of the load on the ATP-bound state is greater than in the ATP-free configuration. This means that 
the force acts as an uncompetitive mixed inhibitor, while in Ref. jl(3] it is shown that kinesin is also 
mixed but competitive. This difference is responsible of the curvature on (v) — F curves at high 
ATP concentration. Figure [2^ shows these curves with the simulation data and the predictions of 
the analytical expression with an excellent agreement. 

Finally, it is interesting to define the coupling ratio and the efficiency and to see how are they 
modified by P. The coupling ratio r can be expressed as the quotient between the total number of 
performed steps and the total number of consumed ATP's. On the other hand, the efficiency 77 can 
be defined as the ratio between the work performed against P, W, and the total input of energy, 
nE. Thus, 

xcm W 
nL nE 
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FIG. 2: a) Simulated mean velocities versus F for p = 1 (circles) and p = 10~ 5 (triangles). The insets show 
V max (F) and Km{F) versus F. Triangles are the values from the table in^ 3 ). Solid lines are the fits of 
Q in order to get Ki U and Ki C . With the two inhibition constants, J5J) can be plotted obtaining the solid 
lines of the main figure, b) Coupling ratio r and efficiency 77. Circles arc simulated data and solid lines are 
predictions of 

where xcm = + ^2) arid supposing that xciu(t = 0) = 0. W = Fxcm, so we can write 
7) = rFL/E. This means that the global efficiency is simply the efficiency in a single step multiplied 
by the coupling ratio. We can go further if we consider Q and the fact that V max is proportional 
to r, and then, 



1 + K iu {F s /F-l) E (1 + Kiu ( Fs -F)> 

FigEb shows the simulated data for r(F) and r/(F) as well as the theoretical predictions. It is 
interesting to remark that the maximum efficiency is slightly below 0.15. 

We have presented an inchworm mechanism which is able to perform directed transport and 
analyzed how it behaves under two variables, the ATP concentration through the probability p 
and the external load F. The motor can be described as a uncompetitive mixed inhibitor obtaining 
analytical expressions for the mean velocity as a function of the two control variables that fits 
accurately the simulated data. Finally, we have discussed the coupling ratio, showing that the 
motor loses the tight coupling as F increases. The efficiency is related with the coupling ratio and 
an analytical expression is given with a good agreement with the simulations. 
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